Effects of running couplings on jet conversion photons 

Lusaka Bhattacharya 

Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Kolkata, Pin-700064, India 



Abstract 

We calculate photons from jet-plasma interaction considering collisional and radiative energy 
loss of jet parton. The phase space distribution of the participating jet is dynamically evolved 
by solving Fokker-Planck equation. We treat the strong coupling constant (a 3 ) as function of 
momentum and temperature while calculating the drag and diffusion coefficients. It is observed 
that the quenching factor is substantially modified as compared to the case when a s is taken 
as constant. It is shown that the Phenix data is reasonably well reproduced when contributions 
from all the relevant sources are taken into account. 
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1. Introduction 

Heavy ion collisions have received significant attention in recent years. Electromagnetic 
probes (photons, dileptons etc) have been proposed to be one of the most promising tools 
to characterize the initial state of the collisions [1] . Because of the very nature of their 
interactions with the constituents of the system they tend to leave the system almost 
unscattered. Photons are produced at various stages of the evolution process. The initial 
hard scatterings (Compton and annihilation) of partons lead to photon production which 
we call hard photons. If quark gluon plasma (QGP) is produced initially, there are QGP- 
photons from thermal Compton plus annihilation processes. Photons are also produced 
from different hadronic reactions from hadronic matter either formed initially (no QGP 
scenario) or realized as a result of a phase transition from QGP. 

These apart, there exits another class of photon emission process via the jet conversion 
mechanism (jet-plasma interaction) [2] which occurs when a high energy jet interacts with 
the medium constituents via annihilation and Compton processes. 

In current heavy ion collision experiments, the temperature T is not only the important 
scale, momentum scale, k, (of the partons) is also important. Therefore running of the 
coupling in the high momentum regime (p ~ T) has to be taken into account to calculate 
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the cross sections and the energy-loss processes. In this work we calculate photons from 
jet-plasma interaction taking into account running of QCD coupling and both collisional 
and radiative energy losses. 

The plan of the article is as follows. We discuss the formalism in the next section. 
Results will be discussed in the section 3. Finally we will conclude. 

2. Formalism 



2.1. Jet- Photon Rate 



The lowest order processes for photon emission from QGP are the Compton scatter- 
ing (q(q) g — > q(q)"f) and annihilation (qq — > 57) process. The differential photon 
production rate for this process is given by [3] : 

E W~p = 2t^Y Jnl1^f J e t (Pi).f2(E 2 )6(p 1 +p 2 -p 3 -p)\M\ 2 (l±ME 3 )) (1) 

where, \M\ 2 represents the spin averaged matrix element squared for one of those pro- 
cesses which contributes in the photon rate and Af is the degeneracy factor of the corre- 
sponding process. fj e t, f 2 and f 3 are the initial state and final state partons. f 2 and f 3 
are the Bose-Einstein or Fermi-Dirac distribution functions. 



2.2. Fokker - Planck Equation: Parton transverse momentum spectra 



In the photon production rate (from jet-plasma interaction) one of the collision partners 
is assumed to be in equilibrium and the other (the jet) is executing random motion in 
the heat bath provided by quarks (anti-quarks) and gluons. Furthermore, the interaction 
of the jet is dominated by small angle scattering. In such scenario the evolution of the jet 
phase space distribution is governed by Fokker-Planck (FP) equation where the collision 
integral is approximated by appropriately defined drag (77) and diffusion coefficients [4]. 

The drag and diffusion coefficients are infrared singular. The infra-red cut-off is fixed 
by plasma effects, where only the medium part is considered, completely neglecting the 
vacuum contribution leading to ambiguity in the energy loss calculation. If the latter part 
is taken into account the strong coupling should be running. Thus for any consistent 
calculation one has to take into consideration this fact. In that case a s = a s (k,T) 
(k = y^\oj 2 — q 2 \ in this case), and the above integrals must be evaluated numerically 
where the infra-red cut-off is fixed by Debye mass to be solved self-consistently: mo{T) = 
4tt (l + ^) a s (mrj(T) 7 T)T 2 Here the strong coupling which we take as running, i. e. 

a s = a s {^\ui 2 — q 2 \, T). We chose the following parametrization of a s which respects 
the perturbative ultra-violet (UV) behavior and the 3D infra-red (IR) point [5]: 

a (k T) = — (2) 

1 + exp(u 2 ± - u 3 ) (l + exp(v 2 ±-v 3 ))(ln(e + (±)« + (±fy {) 



with k = ^/\oj 2 — q 2 in this case. The parameters a, b and A s are given by a — 9.07, 
b = 5.90 and A s = 0.263 GeV. For the limiting behavior (k « T) of the coupling 
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we choose, u\ = a*3d(l + exp(— u 3 )) Here a*^d and a* s denote the values of the IR 
fixed point of SU(3) Yang- Mills theory in d — 3 and d — 4 dimensions, respectively. 
The remaining four parameters (u 2 = 5.47, u 3 = 6.01, v 2 = 10.13 and v 3 — 9.27) fit the 
numerical results for pure Yang-Mills theory obtained from the RG equations in Ref. [6] . 

In our calculation we have considered both collisional and radiative energy losses in 
the following manner. 



V = Vcol\ + ?7rad = ^ 



dE\ fdE 



d x i coil \dx / rad 



(3) 



For running a s , the expressions for the collisional and radiative energy losses can be 
found in [4]. Having known the drag and diffusion, we solve the FP equation using 
Green's function techniques (for details see Ref. [7]). 

2.3. Space time evolution 

In order to obtain the space-time integrated rate we first note that the phase space 
distribution function for the incoming jet in the mid rapidity region is given by (see 
Ref. [8] for details) 

f (~~*t\\ (27r) 3 7?(|uy|) tj l dN 

fjet(r,p,t)\ y=0 = (p T ,t)5(z ) (4) 

VqVU ~ z o Prd 2 p T dy 

With this jet parton phase space distribution function one can easily obtain jet photon 
yield from eqn. (1): 

d/VT f l4 d/VT 
a x 



d 2 pxdy J d 4 xd 2 pTdy 



^LJ C dt'J rdr J dcj>V(w r ) 16(27r ' )7g J dsdt\Mi\ 2 J dE x dE 2 



t, o 



1 dN ^ h(E 2 )(l±f 3 (E 3 )) 

(p 1T ,t)—= = (5) 



Pit dp\ T dy ' y/ a E 2 2 + 2bE 2 + c 

3. Results 



In order to obtain the photon pr distribution we numerically integrate Eq. (5). The 
results for jet-photons for RHIC energies are plotted in Fig. 1 (left) where we have taken 
Ti = 446 MeV and t, = 0.147 fm/c. We find that the yield is decreased with the inclusion 
of both the energy loss mechanisms as compared to the case when only collisional energy 
loss is considered. It is to be noted that when one considers collisional energy loss alone 
the yield with constant a s is more compared to the situation when running a s is taken 
into account (see Fig. 1 left). 

In order to compare our results with high px photon data measured by the PHENIX 
collaboration [9], we have to evaluate the contributions to the photons from other sources, 
that might contribute in this pr range. In Fig. 1 (right) the results for jet-photons corre- 
sponding to the RHIC energies are shown, where we have taken Ti — 446 MeV and U = 
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0.147 fm/c. The individual contributions from hard and bremsstrahhmg processes [10] 
are also shown for comparison. 
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Fig. 1. (color online) Left: pt distribution of photons at RHIC energy. The violet (magenta) curve denotes 
the photon yield from jet-plasma interaction with collisional (collisional + radiative) energy loss. The 
blue curve corresponds to the case without any energy loss and the green curve represents the thermal 
contribution. Right: Our results (the orange line which represents the total photon yield) is compared 
with the Phenix measurements of photon data [9] . 



4. Summary 

We have calculated the transverse momentum distribution of photons from jet plasma 
interaction with running coupling, i. e. with a s = a s (k,T) where we have included both 
collisional and radiative energy losses. Using running coupling we find that the depletion 
in the photon px spectra is by a factor of 2 — 2.5 more as compared to the case with 
constant coupling for RHIC energies Phenix photon data have been contrasted with the 
present calculation and the data seem to have been reproduced well in the low px domain. 
The energy of the jet quark to produce photons in this range (4 < px < 14) is such that 
collisional energy loss plays important role here. It is shown that inclusion of radiative 
energy loss also describes the data reasonable well. 
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